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ABSTRACT 

We present the results of an observational campaign which obtained high time cadence, high precision, 
simultaneous optical and IR photometric observations of three M dwarf flare stars for 47 hours. The 
campaign was designed to characterize the behavior of energetic flare events, which routinely occur on M 
dwarfs, at IR wavelengths to milli-magnitude precision, and quantify to what extent such events might 
influence current and future efforts to detect and characterize extrasolar planets surrounding these stars. 
We detected and characterized four highly energetic optical flares having [/-band total energies of ^7.8 
x 10 30 to ~1.3 x 10 32 ergs, and found no corresponding response in the J, H, or Ks bandpasses at the 
precision of our data. For active dM3e stars, we find that a ~1.3 x 10 32 erg [/-band flare (AU max ~1.5 
mag) will induce <8.3 (J), <8.5 (H), and <11.7 (Ks) milli-mags of a response. A flare of this energy or 
greater should occur less than once per 18 hours. For active dM4.5e stars, we find that a ~5.1 x 10 31 
erg U-band flare (AU ma x ~1.6 mag) will induce <7.8 (J), <8.8 (H), and <5.1 (Ks) milli-mags of a 
response. A flare of this energy or greater should occur less than once per 10 hours. No evidence of stellar 
variability not associated with discrete flare events was observed at the level of ^3.9 milli-mags over 1 
hour time-scales and at the level of ~5.6 milli-mags over 7.5 hour time-scales. We therefore demonstrate 
that most M dwarf stellar activity and flares will not influence IR detection and characterization studies 
of M dwarf exoplanets above the level of ~5-ll milli-mags, depending on the filter and spectral type. 
We speculate that the most energetic megaflares on M dwarfs, which occur at rates of once per month, 
are likely to be easily detected in IR observations with sensitivity of tens of milli-mags. We also discuss 
how recent detections of line flux enhancements during M dwarf flares could influence IR transmission 
spectroscopic observations of M dwarf exoplanets. 

Subject headings: stars: individual (YZ CMi, EV Lac, AD Leo) — stars: flare 



1. INTRODUCTION 

Pioneering exoplanet studies using the excellent sensi- 
tivity and stability of the Spitzer Space Telescope have 
begun to charac terize some of the basic properties of exo- 
planets (see e.g. ISeager fc Deminj [2010 K Secondary tran- 
sit observations have enabled determination of the temper- 
ature s of exoplanets (|Charbonneau et aT]|2005l ; lAgol et ail 
2010), while monitoring of infrared (IR) flux as a func- 
tion of orbital phase, i.e. "phase-mapping", has been 
used to place constraints on the atmospher ic recirculation 
efficiency o f exoplanetary atmospheres (jKnutson et al.l 
120071 I2009D . With respect to exoplanet detections, re- 
evaluations of the impact of stellar activity and flares from 
M dwarfs on the habitability of planets they harbor indi- 
cate that these sys tems could be amenable to hosting hab - 
itable exoplanets (jTarter et al.l l2007t ISegura et al.l 120101) . 
Super-Earth size exoplanets in the habitable zone of M 
dwarfs can be detected via the transit technique using ex- 
isting instrumentation; hence, survey programs such as 
MEarth are targeting large numbers of M dwarfs to find 
these types of planets (|Irwin et al.ll2009D . 

Low-mass stars are well known to exhibit a range of 



stellar activity phenomena (|Hawlevlll993l) . and the preva- 
lence of active s tars is observed to decrease as a function 
of age (see e.g. IWest et al.ll2008l ). One of the more en- 
ergetic byproducts of stellar activity are flares, which are 
thought to be produced by magnetic reconnection in the 
atmospheres of cool stars. Active M dwarfs exhibit some 
of the more energetic and hi g hest rate of flares amongst 
cool stars (lLacv et al.l 119761 : lHawlev fc Pettersenl 1991; 
Kowalski et alll2009l 12010ft . M dwarf flares can last from 
minutes to many hours, exhibit U-band flux enhancements 
of a factor of several hundred above quiescence, and emit 
energies of up to T^u.fiare > 10 34 ergs (lHawley fc Pettersen l 
199ll ICullv et alJll994 iKowalski et alj|2010t lOsten et al.l 
2010). The optical white-light emission has been shown 



to be fit by a two component model comprised of Balmer 
continu um emission plus a T ~10,000K blackbody com- 
ponent (Kowalsk i et al.l [20101 ). Although M dwarf flares 
are most conspi cuous at radio, blue optical, and X - 
ray wavelengths (lHawlev et al.l 120031: lOsten et al.l 120051 ). 
iSchmidt et atl EoTW has recently presented the clear de- 
tection of IR (0.95-2.46^m) flare signatures in H I Paschen 
lines, H I Br7, and the He I 10830 A line during events 
confirmed with simultaneous optical photometry. Several 
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solar flares have also been detected in the continuum re- 
gion near 1.56/xm in the IR, near the H~ opacity min- 
imum, implying formation deep within the solar atmo- 
sphere (|Xu et al.ll200l ICheng et al.ll2010D . 

Some types of stellar activity have already been ob- 
served to "contaminate" exoplanet detection and char- 
acterization studies. Variations in the appearance and 
depth of transit li ght curves have been observed for 
FGKM host stars (TBentlev et all 120091: iDittmann et all 
I2009t lAgol et all l2010t iKundurthv et alj 120111) and Tn 
trans mission spect ra of the M dwarf Gl 4 36 ()Knutson et al.l 
1201 lh . Moreover. iKnutson et all (|2010f ) suggest that UV 
flux from chromospherically active stars may destroy the 
atmospheric compounds responsible for producing temper- 
ature inversions observed in the atmospheres of exoplanets 
orbiting less active stars. 

Although M dwarfs are known to exhibit a range of 
activity, the magnitude of such variability in IR contin- 
uum bandpasses relevant to the detection (via primary 
transits) and characterization (via secondary transits and 
phase mapping) of exoplanets remains poorly understood. 
Marginal variability that may be due to flares was seen in a 
statis tical investigation of sin gle-epoch 2MASS calibration 
data (jDavenport et al.ll201lh . That study is complemen- 
tary to the simultaneous monitoring observations reported 
here, which provide flare rates for individual stars. The 
few previous simultaneous optical-IR observations of M 
dwarf flares re port anomalous IR flux dec rements during 
optic al flares (IRodono fe Cutispotol Il988l ) or null detec- 
tions (jPanagi fe Andrewslll995l ). We caution however that 
the observational techniques and analyses of these simulta- 
neous photometric studies were dubious and not optimized 
for precision IR photometry. 

Could the 5 milli-mag transit signal of a super-Earth 
residing in the habitable zone around a M5V star 
(jlrwin et al.l [2009) be routinely masked by stellar activ- 
ity and flares when observed in the IR? In this paper, we 
present simultaneously obtained, high time cadence opti- 
cal and IR photometric monitoring of well known active 
M dwarfs at milli-mag precision to address this question. 

2. OBSERVATIONS AND DATA REDUCTION 

We used a suite of 3 optical photometric facilities and 1 
IR photometric facility to obtain simultaneous broad-band 
optical and IR photometry of three active M dwarfs known 
to exhibit frequent flares, between 2009 October 2-7 and 
2011 February 11 - 16. The blue optical photometric data 
were obtained to identify the onset and characterize the 
energy of flare events, whereas the red optical photome- 
try was used to help characterize the color of flares. A 
basic summary of the targets observed at each observa- 
tional facility, along with exposure times, filters used, and 
observing cadences is presented in Table [T] Most of our 
2009 observing run was lost due to a tropical storm hitting 
KPNO, while poor weather and instrument failures only 
affected ^2 nights of our 2011 observing run. 

2.1. Optical Photometry 

The Astrophysical Research Consortium Small Aper- 
ture Telescope (ARCSAT) 0.5m at the Apache Point Ob- 
servatory (APO) was one of three facilities we used to 
obtain high candence optical photometry. Observations 



were made with the University of Washington Flare- Cam 
l|Hiltodl20litlHilton et aHl201lD through Sloan g,r,i filters, 
and were recorded using a thermo-electic cooled (TEC) 
1024 x 1024 pixel APOGEE CCD. The full frame readout 
of the 8.' x 8.' FOV was ~1 second. Twilight sky flats and 
dark frames were obtained every night to aid reduction 
of the data. The data were reduced using standard tech- 
niques, and differential aperture photometry was extracted 
using a Python-based script which called the IRAF phot 
package. 

The New Mexico State University (NMSU) automated 
lm telescope at APO also provided high cadence Johnson 
U band observations. Further details a bout the NMSU 
lm an d its operations are described by iHoltzman et al.l 
(|2010h . The data were reduced using an automatic pho- 
tometry pipeline, yielding differential photometry for our 
science targets. 

During our 2009 observing run we also obtained simul- 
taneous optical monitoring of our science targets with the 
Wisconsin Indiana Yale NOAO (WIYN) 0.9m telescope. 
We used the S2KB CCD, a 2048 x 2048 pixel camera with 
a pixel scale of 0'.'60 pixel -1 . We employed 2x2 bin- 
ning of the S2KB chip, and a 260 x 260 pixel sub-array of 
the full frame, to reduce read out time and obtain higher 
cadence photometry. The data we present were obtained 
in the Harris U band filter. Nightly bias and dome flat 
field frames were obtained, and the data were reduced us- 
ing standard IRAF techniques. Differential photometry 
was extracted using the same Python-based software im- 
plemented for our ARCSAT data. 

2.2. IR Photometry 

Our IR photometry was obtained at the KPNO 2.1m 
telescope using the Simultaneous Quad Infrared Imaging 
Device (SQIID). SQIID obtains simultaneous J (X cen ter 
= 1.267 /mi), H {Xcenter = 1-672 (im), and Ks (A centej . = 
2.224 fj,m) band imagery in 3 of the 4 quadrants of its ar- 
ray. At the f/15 focus, each 512 x 512 quadrant of SQIID's 
ALADDIN array spans a field of view of 5'07 x 5.'28 with 
a pixel scale of 0.69 arcsec pixel -1 . Our science exposures 
were generally obtained at or near the array's minimum 
exposure time of 0.874 seconds, and we defocussed the 
telescope until stars took the appearance of ~20? - ~28'.' 
-wide donuts, to keep the counts below the array's linear- 
ity limit. Because of large overheads in the transmission 
of data from the array, our net observing cadence was ~60 
seconds (Table [T]). 

Night-time, deep blank-sky images were obtained every 
clear night to serve as sky flats; dark frames were also ob- 
tained every day corresponding to the integration times 
of our sky flats. We used a 3-4 dither point pattern in 
order to mitigate ghosts produced from observing bright 
sources. Since we were interested in performing differen- 
tial photometry, the specific number and pattern of dither 
points we adopted varied by source, to ensure each frame 
contained our science target and at least 1 comparison 
star. We periodically adjusted the position of the science 
target on the array to ensure that small amplitude, long 
time-scale telescope tracking imperfections did not cause 
the target to fall on different pixels, thus making it more 
susceptible to imperfections in the flat fielding process. 

We developed a custom data reduction routine to ex- 
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tract the highest precision photometry possible. We first 
applied a third order linearity correction to the data using 
coefficients in the SQIID user's manual, and then divided 
these data by a normalized flat field image. Sky subtrac- 
tion and dark correction of our science images was achieved 
by subtracting immediately adjacent dither frames from 
each other. Individual pixels >16,000 ADUs (the linearity 
limit) or < -500 ADUs were next flagged as bad or warm 
pixels, and corrected using the IRAF routine fixpix. Next, 
low-level array artifacts, such as column banding and less 
frequently observed row banding, were characterized by 
sampling this structure in star-free regions of the array 
and removed via custom IDL routines. 

To obtain photometry for our SQIID data, we first 
identified the centr oid position of each star in our FOVs 
using SExtractor (|Bertin fc Arnoutsl [l996), and a cus- 
tom reference profile which mimicked the extended donut- 
shape of our sources. Next, we ran these stellar posi- 
tions through IRAF's phot routine to extract aperture 
photometry, using a target aperture of radius 20 pixels 
(EV Lac and YZ CM) to 24 pixels (AD Leo) and a 4 
pixel wide background sky annulus starting at a radius of 
30 pixels from each centroid position. We then extracted 
differential photometry for our science targets using the 
brightest, and often only, comparison star with the SQIID 
FOV, as summarized in Table [2] Previous IR monitor- 
ing progr ams dedicated to o bserving M dwarfs for transits 
(MEarth; E rwm eTal1l20ll or high precision observations 
of brown dwarfs to character ize photospheric condensates 
(Bail er- Jones fe Lamml [2~003f) have described how varia- 
tions in relative humidity can affect IR differential pho- 
tometry. We found we were able to remove most of the 
small amplitude, long time-scale photometric trends in our 
data likely attributable to these effects by simply fitting 
a linear function to each dither position's light curve in- 
dividually. Finally, each dither position was normalized 
to a uniform scale and combined to yield the highest time 
cadence differential photometry from our observations. 

Table [3] summarizes the relative photometric stability 
we were able to obtain for every target and filter, both 
across an entire night and across 1-hour windows, quanti- 
fied as the observed standard deviation. The IR stability 
we achieved for AD Leo, 9.7-12.6 milli-mags over a full 
night and 8.4-10.5 milli-mags over a 1 hour window, was 
systematically larger than that achieved for our other tar- 
gets due to the lower flux of the brightest comparison star 
in its field of view. The best IR photometric stability we 
achieved was an impressive 5.1-6.0 milli-mags over a full 
night and 3.8-3.9 milli-mags over a 1 hour window for our 
2011 February 13 monitoring of YZ CMi. 

3. RESULTS 

Figure Q] shows a representative multi-filter differential 
photometric light curve obtained for the star AD Leo on 
2011 February 13. Note that photon statistics-based er- 
ror bars are plotted in all panels of Figure [TJ and for all 
figures in this paper, but these error bars are generally 
smaller than the size of the data points. Before charac- 
terizing the behavior of stellar flares in our IR data, we 
first identified these events in our simultaneously obtained 
optical photometry. We used t h e IDL -bas ed flare find- 
ing software described in iHiltonl (|201ll) and iHilton et al.l 



(|201lf) . Flares are identified as a single epoch which is 3.5- 
(7 brighter than the standard deviation of the local mean, 
followed by three epochs which are 2-cr above the local 
mean. Each flare identified was also double-checked by 
eye to confirm the event was robust. The total number 
of flares detected in each bandpass using this flare finding 
algorithm is shown in Table [3] a total of 20 flares were 
detected in our U -band data. 

For the larger optical flare events observed, additional 
flare properties such as the peak magnitude enhancement 
and total flare duration are compiled in Table |3J Flare 
equivalent durations, defined as the amount of time each 
star would need to spend at a quiescent level to produce 
the s ame total energy as during each flare (|Gershberel 
|1972|) . were determined via our flare finding software by 
integration under each light curve. We adopted the qui- 
escent luminosities for AD Leo (1.59 x 10 29 ergs s _1 ), EV 
Lac (6.63 x 10 28 e rgs s" 1 ) and YZ C Mi (4.57 x 10 28 erg s 
s" 1 ) computed by IHiltonl (1201 lD a ndlHilton et al.l (<201l[ i. 
and based on data in iReid et all (119951) . We then multi- 
plied the quiescent luminosities by the equivalent duration 
of each event, to extract flare energies (Table 01 . 

3.1. Optical Flares 

As the focus of this paper is on characterizing the prop- 
erties of stellar flares in the IR, we focus our attention 
on the four largest U-band flares we observed on AD Leo 
(event #1), YZ CMi (events #2 & #3), and EV Lac 
(event #4) , as listed in Table 0] and illustrated in Fig- 
ure [3J Our flares exhibited U -band energies ranging from 
~7.8 x 10 30 (event #4) to -1.3 x 10 32 (event #1) ergs . 
Flare events #1-3 are considered strong (jLacv et al. ll976) 
although the y are less energe t ic tha n the very rare ~10 34 

(f2010h obs e rved o n YZ CMi. 



erg megaflare lKowalski et al 



For each of our targets, 



Lacv et all (|1976l) presented 



flare frequency distributions which follow power la ws, ex- 
cept a t th e lowest and highe st energy regimes. IHiltonl 
(|201ll) and lHilton et al.l (|201ll ) expanded upon this earlier 
work, presenting flare frequency distributions for a wider 
range of M dw arfs. After correcting for detection efficiency, 
Hilton (2011) was able to demonstrate that flare frequency 
distributions follow power law distributions even down to 
low ener gies. W e use the fl are frequency distributions pre- 
sented in lHiltonl (|201ll) and lHilton et al.l (|201lD for this pa- 
per. Note that since the flare energies we computed in Ta- 
ble HI used the same quiescent stellar luminosities employed 
to construct the Hilton et al flare frequency distributions, 
any future refinements of these quiescent luminosities will 
not influence the flare rates we hereafter cite. 

We determine that a flare having equal or greater U- 
band energy to our event #1, on AD Leo, would occur 
less than once per 18 hours. For YZ CMi, we find that 
a flare having equal or greater U -band energy for events 
#2-3 would occur less than once per 10 hours, and the 
EV Lac flare frequency distribution indicates that a flare 
having equal or greater £/-band energy to our event #4 
would occur less than once per 7 hours. 

3.2. IR Continuum Flares 

We utilized the same flare finding algorithm described 
in Section [3] and used to identify flares in our optical data 
to search for contemporaneous flares in our full set of IR 
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data. We detected no statistically significant deviations in 
our IR photometry at the epochs of the four major, op- 
tically detected flares. We quantify the 1-er upper limits 
of these four events by computing the standard deviation 
of all data in a 40 minute time window before each event, 
and summarize these limits in Table 2] 

Our [/-band photometry provides us with information 
about the starting time of each flare event. As seen 
in Figure [J and demonstrated by many othe rs (see e.g. 
lHawlev fc Pettersenl Il99lt Iffilton et alll2011l) . the total 
duration of flares in the optical is color dependent. Flare 
#1 in our survey, for example, was visible in the U filter for 
lOx longer than it was in the i filter (Table [4j. We there- 
fore estimated upper limits to the duration which each 
flare event would be visible in the J, H, and Ks filters by 
adopting the duration of each event in the reddest optical 
filter in which it was detected (e.g. the i-band for event 
#1). Our observing cadence (Table [T]) constrains the un- 
certainty in the durations we quote. This information en- 
abled us to time bin our data, as shown in Figures|«2S and 
[5]for the J, H, and Ks filters respectively, to intervals which 
would create two data points (i.e. ~2 minute binning) 
and 1 data point (i.e. ^4 minute binning) across the up- 
per limit duration of each event. As flare light curves can 
exhibit significant changes over 2-4 minute intervals (see 
e.g. Figure [1]), binning our data to these longer cadences 
could cause flare signatures to be diminished. Conversely, 
binning will help to reduce the dispersion of quiescent- 
phase data, and therefore could help elucidate the pres- 
ence of low-amplitude slowl y varying even t s, suc h as the 
long-duration flare shape of lHawlev et al.l (|l995l ). How- 
ever, even the binned data did not reveal a statistically 
significant event in the IR. 

Following practices described in iKowalski et al.l (|2009f ), 
we also computed a We lch and Stetson variability index 
( Wel ch fc Stetson! Il993l ). $jhk, for each epoch of our 
data. This type of variability index builds on the prin- 
ciple that the photometric variations we are searching for 
should appear at similar time epochs across multiple fil- 
ters; hence, it serves as a way to identify faint signals com- 
mon to each of our J, H, and Ks filter data. In practice, 
the index is computed simply by: 

/ Flux j\ ( Flux H \ ( Flux Ks \ „> 

" JHK = l^rj {—) l^rj (1) 

where the flux in each multiplicative term is deter- 
mined from each filter's differential photometry light 
curve. The resultant <&jhk index for our unbinned, ^2 
minute binned, and ~4 minute binned data exhibited no 
statistically significant change during the time of each flare 
event, as shown in Figure |5] for event #1, as compared to 
the value of the index preceeding and following each event. 

We can not exclude the possibility that our IR observing 
cadence (Table [1} of one (0.874 second) integration every 
~60 seconds does not influence our ability to detect the 
optically observed flare events. To help better quantify 
whether our observations could miss IR responses to flares 
which begin and end on time-scales faster than the relative 
observing cadence we achieved, we identified the time of 
the peak flare emission during each event, as listed in Table 
[S] The quoted error bars on these times correspond to the 
epoch of the prior and subsequent integration in each filter. 



We also list the epochs of the nearest 1-2 IR observations 
in Table [5] relative to the observed epoch of the peak AU- 
band flux. The epochs of our IR observations essentially 
overlap with the epochs of the observed peak A [/-band 
flux for flare events #1,3, and 4 to within the timing un- 
certainties. Our closest IR integration for flare event #2 
was 18.5 seconds after the oberved AU-band peak, which 
is larger than the ±13 second uncertainty in the epoch of 
the peak A U -band flux. These data quantify the limits at 
which our different optical versus IR observing cadences 
could have resulted in the non-detection of flare events. 

We conclude that, at the limit of our observational data, 
we find no compelling evidence of flare enhancements in 
our J, H, or Ks filter data during the ~7.8 x 10 30 to ~1.3 
x 10 32 erg [/-band flares we detected. 

4. DISCUSSION 

We have presented the results of ~47 hours of high ca- 
dence, high precision, simultaneous optical and IR pho- 
tometric monitoring of 3 active M dwarfs. We now dis- 
cuss the interpretation of these results in the context of 
the effects that stellar flares could have on a variety of 
M dwarf exoplanet studies. For reference, we remind the 
reader that a Super-Earth (2 R^) residing in the habitable 
zone around a dM5 star (~0.07 4 AU separation; see e.g. 
iNutzman fc Charbonneaull2008l ) would have an orbital pe- 
riod of 14.5 days and a transit duration of ~2 hours. 

4.1. Effects of Stellar Flares on IR Transit and Phase 
Mapping of M Dwarf Exoplanets: Observational 
Constraints 

In Section [32] we demonstrated that during four signif- 
icant optical flare events, having U -band energies ranging 
from ^7.8 x 10 30 to ~1.3 x 10 32 ergs, we observed no sta- 
tistically significant evidence of corresponding broad-band 
enhancements in the J, H, and Ks filters at the 5.1-11.7 
milli-mag level. These upper limits were computed from 
the standard deviation of all data in a 40 minute time 
window before each flare. 

We quantified the relative frequency that one would ex- 
pect for flares with energies similar to events #1-4 in Sec- 
tion 13.11 using flare freq uency d i stribu tions observed for 
our sp ecific dMe stars in iHiltonl (|2011h and iHilton et al.l 
(|2011|) . By combining these two observational properties, 
we can place upper limits on the effects of stellar flares 
on future continuum-based observations of M dwarf exo- 
planetary systems. For active M3Ve stars, we find that a 
~1.3 x 10 32 erg [/-band flare will induce <8.3, <8.5, and 
<11.7 milli-mags of an effect in the J, H, and Ks filters 
respectively. A flare of this energy or greater should occur 
less than once per 18 hours. For active M4.5e stars, we 
find that a ~5.1 x 10 31 erg [/-band flare will induce <7.8, 
<8.8, and <5.1 milli-mags of an effect in the J, H, and Ks 
filters respectively. A flare of this energy or greater should 
occur less than once per 10 hours. Moreover, we observe 
no evidence of stellar variability not associated with dis- 
crete flare events at the level of >3.9 (J), >3.8 (H), and 
>3.9 (Ks) milli-mags over 1 hour time-scales and no level 
of stellar variability at the level of>6.0 (J), >5.6 (H), and 
>5.1 milli-mags over 7.5 hour time-scales. These data pro- 
vide quantitative upper limits to the level of broad-band 
stellar variability which could be expected in IR transit 
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and phase mapping observations of M dwarfs. 

4.2. Effects of Stellar Flares on IR Transit and Phase 
Mapping of M Dwarf Exoplanets: Theoretical 
Constraints 

Although multi-wavelength radiative hydrodynamic 
modeling of solar and d Me flares has been constructed 
(|Allred et al.ll2005ll2006l ). we are not aware of any suite of 
detailed low mass star flare mode ls in the existing litera- 
ture that present model IR fluxes. iDavenport et al. (|201lh 
begins to mitigate this deficienc y by adopting the tw o 
component flare emission model of Kowalsk i et al.l (|201df ) , 
comprised of Balmer continuu m predictions by radiative 
hydrodynamic flare models of ()Allred et al.ll200a) and a T 
~10,000K blackbody, and extrapolating this model to pre- 
dict expected fluxes over a range of optical and IR filters. 
For our flare event #1 (Table[4}, this model predicts a peak 
j-band flux enhancement (^0.02 mag) which is within a 
factor of 2 of our observed value of 0.04 magnitudes and 
predicts an IR response of 3.5 milli-mags (J-band), 2.1 
milli-mags (if-band), and 1.8 milli-mags (ifs-band) which 
is consistent with the upper limits we set observationally 
in Section [3.21 

Using the IDavenport et al.l (|2011[) extrapolations, we 
explore the level of flare which would need to occur to 
be detected by our observing program, given the typi- 
cal standard deviation of our IR data during epochs of 
observed flares, ~8 milli-mags. For a dM3e star, these 
models suggest a flare exhibiting a peak U -band enhance- 
ment of 2.6 magnitudes (i.e. one magnitude greater than 
our largest flare) would have produced J-band (11.6 milli- 
mags), if -band (6.8 milli-mags), and ifs-band (5.7 milli- 
mags) responses which would have been detected by our 
observational program. Constraining the frequency of such 
large flares is difficult as observed flare frequency dis- 
tribu tions deviate from simple power laws at high ener - 
gies (|Lacv et al.1 fl976l: Iffiltonl [20lTL iHilton et al.l l20"ll . 
Nonetheless, we can conservatively state that the fre- 
quency of such events should be less than once per 35 
hours. 

We also note that the Ef/ >10 34 er gs meg aflare 
on YZ CMi reported by IKowalski et al.l (|2010h was 
also detected by the MEarth monitoring program 
(Irw in et al. 201 If), which obser ves with an i+z filter 
dNutzman fc Charbonneaul l2008f ). This ~6 magnitude 



U band event (Kowalski et al 



2010) produced a -0.6 



magnitude enhancement in the MEarth (i+z) filter (Ir- 
win 2011; personal communication), which is generally 
consistent with the pre dicted level of enhancement from 
IDavenport et al.l (|20 1 it) . For that flare, the model pre- 
dicts J-band (119 milli-mag), £f-band (77 milli-mag), and 
Ks-h&nd (62 milli-mag) enhancements, which are signifi- 
cantly larger than our photometric noise floor. Although 
we did not observe such a large flare, which IKowalski et al.l 
( 2010) suggest should occur at a rate of once per month, 
we speculate that such events will easily be detectable in 
future IR observations of M dwarfs. 

The IKowalski et al.l (|2010f) observational and modeling 
work whi ch formed the basis of the extrapolations com- 
puted by IDavenport et all (|2011[ ) is unique in that it pro- 
vided convincing evidence that the blue optical proper- 
ties of flares could be reproduced by 2 components, a 



Balmer conti nuum predicted b y radiative hydrodynamic 
flare models (jAllred et al.ll2006D and a T ~10,000K black- 
body. This framework represents an important step to- 
wards advancing our understanding of flares; however, it 
is phcnomenological in nature and by design doesn't re- 
produce the entire myriad of obs ervational characteris- 
tics of flares. Kowals ki et al.l (|201ll ). for example, suggest 
that additional components from H I Paschen continuum 
emission and photospheric backwarming can reproduce the 
residual observed flux from 5000-5500 A, which isn't ac- 
counted for in the 2-component computation. It is there- 
fore plausible that including additional components rele- 
vant to the J,H,Ks wavelength regime, such as Brackett 
and Pfund continuum emission, could also be important. 
Although extrapolating a blue/near- UV continuum model 
may not yield a complete description of the IR flare con- 
tinuum, we sugg e st ext rapolations like those presented in 
IDavenport et al. (|201l[) and adopted here can provide an 
order of magnitude guidance on the expected IR response 
of flares. Although beyond the scope of this paper, future 
models of stellar flares should consider providing broad- 
band flux predictions at IR wavelengths to the community, 
to complement the observational constraints provided by 
this work. 

4.3. Effects of Stellar Flares on Transmission 
Spectroscopy Studies of M Dwarf Exoplanets: 
Observational Constraints 

Since the discovery of sodium in th e atmosphere of HD 
209458b bv lCharbonneau et alj (|2002l) . transmission spec- 
troscopy has been recognized as one of the best ways to 
directly probe the chemistry of exoplanetary atmospheres. 
The technique has been successfully employed with cur- 
rent ground- and space-based facilities, leading to detec- 
tions of atmospheric molecules like water, methane, ca r- 
bon monoxide, and carbon dioxide (|Swain et al.ll2009alrb1 . 
l2010t iKnutson et al.l 1201 lh . and is anticipated to play an 
important role in characterizing the atmospher es of ex- 
oplanets once JWST is operational (see e.g. IClampirj 
2010 ; iBelu et al.l 1 2 1 lh . In this context, we note that 
Schmidt et al.l (| 2 1 lh reported the first detection of IR line 
emission in H I Br7 during a ^4 x 10 32 erg (it-band energy) 
flare on EV Lac, whereby the line flux was e nhanced during 
the fla re by ~5% above quiescent values. ISchmidt et al.l 
(|201lh also reported the presence of H I Pa/3, Pa7, and 
Pa<5 emission during three separate flare events with u- 
band (and also [/-band) energies ~4 x 10 31 to ~4 x 10 32 
ergs, with observed line flux enhancements of 5-20% above 
that observed in quiescence. The expected rate, i.e. duty 
cycle, of IR line emiss ion was estimated to be ~3% for ac- 
tive mid-M dwarfs bv lSchmidt et al.l (|201lh . The widths of 
these lines are much smaller than the widths of molecular 
features being characterized in current exoplanet studies, 
and therefore should not preclude investigations of the at- 
mospheric molecular content of exoplanets surrounding M 
dwarfs. However, efforts to interpret future transmission 
spectroscopy observations of M dwarfs should allow for 
the possibility of variable stellar H I line flux at the 5-20% 
level, with a duty cycle of ~3%. If observed in broad- 
band filters, these line flux enhancements would produce 
enhancements of 0.4 milli-mags (J-band) to 0.3 milli-mags 
(iTs-band), significantly below the upper limits of the IR 
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photometric stability presented in Section [3.21 

5. CONCLUSIONS 

We have presented ~47 hours of high cadence, high pre- 
cision, simultaneous optical and IR photometric monitor- 
ing of 3 active M dwarfs. We detected and characterized 
four highly energetic optical flares having U-band total 
energies of ^7.8 x 10 30 to ~1.3 x 10 32 ergs, and found no 
corresponding response in the J, H, or Ks bandpasses at 
the precision of our data. To summarize our results: 

• For active dM3e stars, we find that a ^1.3 x 10 32 
erg [/-band flare will induce <8.3 (J), <8.5 (H), 
and <11.7 (Ks) milli-mags of a response. A flare of 
this energy or greater should occur less than once 
per 18 hours. 

• For active dM4.5e stars, we find that a ~5.1 x 10 31 
erg [/-band flare will induce <7.8 (J), <8.8 (H), 
and <5.1 (Ks) milli-mags of a response. A flare of 
this energy or greater should occur less than once 
per 10 hours. 

• No evidence of stellar variability not associated with 
discrete flare events was observed at the level of ^3.9 
milli-mags over 1 hour time-scales and at the level of 
~5.6 milli-mags over 7.5 hour time-scales. We there- 
fore demonstrate that most M dwarf stellar activity 
and flares will not influence IR detection and charac- 
terization studies of M dwarf exoplanets above the 
level of ~5-ll milli-mags, depending on the filter 
and spectral type. 

• Based on the detection of the Ejj >10 34 ergs 
megaflare event (|Kowalski et al.lfeOlOl ) in the (i+z) 



filter of the MEarth survey (jlrwin et al.ll201ll) . and 
extrapolations of phenomenological flare models to 
IR (J,H,Ks) wavelengths, we speculate that the 
most energetic megaflares on M dwarfs, which occur 
at rates of once per month, are likely to be easily de- 
tected in IR observations which obtain sensitivities 
of tens of milli-mags. 

• Future IR transmission spectroscopic studies of M 
dwarf exoplanets should note that 5-20% flux en- 
hancements have been observed in H I Paschen and 
Br7 lines during flares having total energies similar 
to event #1 in our survey. This level of line flux will 
only produce a small (0.3-0.4 milli-mag) enhance- 
ment in broad-band (J,H,Ks) filter observations and 
should not preclude molecular line analysis. 
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Table 1 
Summary of Observations 



UT Date 




Science Object 


Facility 


Filter 


Integration Time 

sec 


Total Obs Time 
hours 






Obs Cadence 1 






2009 


Oct 


8 


EV Lac 


WIYN 0.9m 


U 


15 


8 






1 obs per 51 sec 






2009 


Oct 


8 


EV Lac 


KPNO 2.1m 


J,H,Ks 


0.87 


8 






1 obs per 52 sec 






2011 


Feb 


12 


YZ CMi 


NMSU lm 


U 


10 


8.25 






1 obs per 18 sec 2 






2011 


Feb 


12 


YZ CMi 


ARCSAT 0.5m 


g 


1 


8.25 


2 


obs 


over 4.2 sec; then 17 


sec 


break 


2011 


Feb 


12 


YZ CMi 


ARCSAT 0.5m 


r 


0.5 


8.25 


3 


obs 


over 4.8 sec; then 16 


sec 


break 


2011 


Feb 


12 


YZ CMi 


ARCSAT 0.5m 


i 


0.5 


8.25 


4 


obs 


over 6.0 sec; then 15 


sec 


break 


2011 


Feb 


12 


YZ CMi 


KPNO 2.1m 


J,H,Ks 


0.87 


8.25 






1 obs per 60 sec 






2011 


Feb 


12 


AD Leo 


NMSU lm 


U 


4 


3 






1 obs per 12 sec 2 






2011 


Feb 


12 


AD Leo 


ARCSAT 0.5m 


g 


1 


3 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


12 


AD Leo 


ARCSAT 0.5m 


r 


1 


3 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


12 


AD Leo 


ARCSAT 0.5m 


i 


1 


3 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


12 


AD Leo 


KPNO 2.1m 


J,H,Ks 


0.87 


3 






1 obs per 58 sec 






2011 


Feb 


13 


YZ CMi 


NMSU lm 


U 


10 


7.5 






1 obs per 18 sec 2 






2011 


Feb 


13 


YZ CMi 


ARCSAT 0.5m 


g 


1 


7.5 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


13 


YZ CMi 


ARCSAT 0.5m 


r 


1 


7.5 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


13 


YZ CMi 


ARCSAT 0.5m 


i 


1 


7.5 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


13 


YZ CMi 


KPNO 2.1m 


J,H,Ks 


0.87 


7.5 






1 obs per 60 sec 






2011 


Feb 


13 


AD Leo 


NMSU lm 


U 


4 


3 






1 obs per 12 sec 2 






2011 


Feb 


13 


AD Leo 


ARCSAT 0.5m 


g 


1 


3 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


13 


AD Leo 


ARCSAT 0.5m 


r 


1 


3 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


13 


AD Leo 


ARCSAT 0.5m 


i 


1 


3 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


13 


AD Leo 


KPNO 2.1m 


J,H,Ks 


0.87 


3 






1 obs per 58 sec 






2011 


Feb 


15 


YZ CMi 


NMSU lm 


U 


10 


9 






1 obs per 18 sec 2 






2011 


Feb 


15 


YZ CMi 


ARCSAT 0.5m 


g 


1 


9 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


15 


YZ CMi 


ARCSAT 0.5m 


r 


1 


9 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


15 


YZ CMi 


ARCSAT 0.5m 


i 


1 


9 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


15 


YZ CMi 


KPNO 2.1m 


J,H,Ks 


0.87 - 1.200 


9 






1 obs per 58 sec 






2011 


Feb 


16 


YZ CMi 


NMSU lm 


U 


10 


8 






1 obs per 18 sec 2 






2011 


Feb 


16 


YZ CMi 


ARCSAT 0.5m 


g 


1 


8 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


16 


YZ CMi 


ARCSAT 0.5m 


r 


1 


8 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


16 


YZ CMi 


ARCSAT 0.5m 


i 


1 


8 


4 


obs 


over 8.4 sec; then 27 


sec 


break 


2011 


Feb 


16 


YZ CMi 


KPNO 2.1m 


J,H,Ks 


0.87 


8 






1 obs per 58 sec 







Note. — The integration times and effective observing cadence are tabulated for each of our observations. Note that the observing cadence 
includes exposure times as well as read-out and data transfer overheads. 1 All cadences compiled in column 7 are approximate. 2 Indicates 
the most common observing cadence, although the robotic NMSU lm telescope occasionally intersperses a single exposure separated by ~24 
seconds, and routinely pauses to refocus roughly every 60 exposures. 
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Table 2 

Summary of Science and Comparison Stars 



Science Object 


Spectral Type 


H Magnitude 


IR Comparison Star 


Spectral Type 


H magnitude 


YZ CMi 


c!M4.5e 


6.005 


HD 62525 


G5 


6.035 


AD Leo 


dM3e 


4.843 


TYC 1423-165-1 


G5 


7.686 


EV Lac 


dM3.5e 


5.554 


BD+43 4303 


KOV 


6.685 



Note. — Basic properties of our three science targets and the associated comparison stars used for computing differential photometry for 
each is compiled. 
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Table 3 
Summary of Detected Flares 



UT Date 


Star 


Filter 


# Flares Detected 


IR Stability (full night) 


IR Stability (1 hr) 










milli-mag 


milli-mag 


20091008 


EV Lac 


U 


1 










.1 





9.0 


5.7 






H 





8.9 


4.5 






Ks 





8.7 


6.5 


20110212 


YZ CMi 


U 


6 










g 


1 










r 


1 










i 













J 





8.6 


6.0 






H 





8.7 


6.1 






Ks 





7.4 


6.1 


20110212 


AD Leo 


U 













g 













r 













i 













.1 





10.0 


9.6 






H 





9.7 


8.0 






Ks 





12.6 


10.5 


20110213 


YZ CMi 


U 


1 










g 













r 













i 













.1 





6.0 


3.9 






H 





5.6 


3.8 






Ks 





5.1 


3.9 


20110213 


AD Leo 


U 


2 










g 


2 










r 


1 










i 


1 










.1 





11.4 


8.5 






H 





11.2 


8.4 






Ks 





12.1 


8.8 


20110215 


YZ CMi 


U 


1 










g 













r 













i 













J 





7.8 


4.5 






H 





8.8 


5.0 






Ks 





8.3 


4.5 


20110216 


YZ CMi 


U 


3 










g 


2 










r 













i 













.1 





6.1 


4.2 






H 





6.2 


4.8 






Ks 





8.9 


5.0 



Note. — The number of flares detected via the automated flare finding algorithm discussed in Section \3\ is summarized for each filter. For 
the SQIID IR (JHKs) photometry, we also compile the observed standard deviation of the differential photometry across the entire night and 
in a 1-hour subset of each night, in units of milli-magnitudes. 
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Table 4 
Summary of Flare Properties 



UT Date 


Star 


Event # 


Filter 


Peak A Magnitude 


Duration 


Equiv Duration 


Flare Energy 


Flare Frequency 










magnitude 


sec. 


sec. 


ergs 


hours _1 


20110213 


AD Leo 


1 


U 


1.53 


2460 


797 


~1. 27c+32 


~0.056 








g 


0.25 


1085 


54 












r 


0.10 


921 


21 












i 


0.04 


232 


4 












J 


< 0.0083 
















11 


< 0.0085 
















Ks 


< 0.0117 










20110212 


YZ CMi 


2 


U 


1.61 


3574 


1121 


~5.12c+31 


-0.100 








g 


0.30 


621 


52 












r 


0.09 


445 


17 












i 


< 0.0227 
















J 


< 0.0078 
















11 


< 0.0088 
















Ks 


< 0.0051 










20110213 


YZ CMi 


3 


U 


1.03 


3973 


722 


~3.30e+31 


-0.128 








.1 


< 0.0067 
















H 


< 0.0070 
















Ks 


< 0.0057 










20091008 


EV Lac 


4 


U 


0.31 


901 


118 


~7.82e+30 


-0.142 








J 


< 0.0072 
















H 


< 0.0102 
















Ks 


< 0.0085 











Note. — The peak delta magnitude, duration, equivalent duration, energy, and flare frequency is given as a function of filter for each of 
our four flare events. Note that positive peak delta magnitudes correspond to flux enhancements. Section 13. 1 1 describes how flare equivalent 
durations, energies, and rates were computed. All quoted peak magnitude upper limits represent the 1-<t standard deviation of the observation 
data. 
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Table 5 

Timing of Optical and IR Flares 



UT Date 


Star 


Event # 


Filter 




Relative Time from Flare Peak 


20110213 


AD Leo 


1 


U 




0±10.4 sec 








g 




21.5±1.6 sec after U-filter 








r 




2.4±1.6 sec after U-filter 








i 




9.7l2 5 6 5 sec after U-filtcr 








J,H,Ks 


11 


sec before U-filter; 46 sec after U-filtcr 1 


20110212 


YZ CMi 


2 


U 




0±13 sec 








g 




3.7tjg 6 8 sec after U-filter 








r 




1.5 t, 5 fi 8 sec after U-filter 








J,H,Ks 


42 sec before U-filtcr; 18.5 sec after U-filtcr 1 


20110213 


YZ CMi 


3 


U 




0±13 sec 








J,H,Ks 


12 


sec before U-filter; 48 sec after U-filtcr 1 


20091008 


EV Lac 


4 


U 




± 43.5 sec 








J,H,Ks 




1.4 sec after U-filter 1 



Note. — We identified the relative epoch of peak flare emission in every filter that each of our four flare events were detected. The quoted 
error bars on these time correspond to the time of the prior and subsequent integration. x We did not detect clear evidence of flares in the 
J,H,Ks filters; thus, for these filters we list the nearest 1-2 epochs of exposures relative to epoch of the peak A magnitude flux observed in the 
U-filter. 
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AD Leo - 20110213 



< 



-1.5 
-1.0 
-0.5 
0.0 
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-0.3 
-0.2 
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-0.2 
u -0.1 
< 0.0 
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-0.05 
3 0.00 
0.05 



-0.05 =- 
0.00 
0.05 



< 



< 



-0.05 
0.00 
0.05 

-0.05 
0.00 
0.05 



••••• .• 



.--•\ 



10.0 



10.5 



11.0 



11.5 



12.0 



12.5 



UT (hours) 



Fig. 1. — The differential photometric light curve for our observations of AD Leo on the night of 2011 February 13 is shown for the U, g, r, 
i, J, H, and Ks filters. Photon statistics-based error bars are plotted in all panels, but these error bars are generally smaller than the size of 
the data points 
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Fig. 2. — A zoomed view of the large flare observed on AD Leo on 2011 February 13 (Figure[T}, referred to as flare event #1 in Table|4] 
The online version of this journal contains analogous figures for flare events #2 - #4. 
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Fig. 3. — The J-band differential photometry of flare event #1 (Table |4). observed on AD Leo on 2011 February 13. The ~4 minute 
duration of the flare in the i filter, the closest filter to the J-band, is depicted by solid vertical lines. The top panel depicts the native time 
resolution of the data (~1 minute), the middle panel depicts the data binned to a cadence of ~2 minutes, and the bottom panel depicts 
the data binned to a cadence of ~4 minutes. The online version of this journal contains analogous J-band figures for flare events #2 - #4 
described in Table [4] 
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Fig. 4. — The same as Figure [3] plotted instead in the H-band. The online version of this journal contains analogous H-band figures for 
flare events #2 - #4 described in Table [4] 
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Fig. 5. — The same as Figure[3] plotted instead in the Ks-band. The online version of this journal contains analogous Ks-band figures for 
flare events #2 - #4 described in Table [4] 



17 



Event 1 



2.0xlO e 



1.6x10° ~. 



1.2xlO e 















• • •• • 


• ■ 






• • .- • • • .*• 




— • 

















10.0 



10.5 



11.0 



11.5 



12.0 



12.5 



-2 minute bin 



3.8x10° E- 



ae 



3.4x10° ^. 



3.0x10 



6 



10.0 



10.5 



11.0 



11.5 



12.0 



12.5 



ac 



1.2x10 
1.1x10 

l.oxio'l ' 

9.0x10 



10.0 



10.5 



~A minute bin 




11.0 11.5 
UT (hours) 



12.0 



12.5 



Fig. 6. — The Qjhk statistic for flare event #1, depicted by vertical lines, is shown. We find no statistical evidence for an enhanced signal 
in this statistic for our unbinncd data (top panel), our data binned in ~2 minute intervals (middle panel), or in our data binned in ~4 minute 
intervals (bottom panel). 



